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Due to the lack of both parity and time reversal symmetries, the Rashba superconductors CePtsSi, 
' CeRhSi3, and CeIrSi3, in the presence of a magnetic field, are unstable to helical (single plane- 

wave) order. We develop a microscopic theory for such superconductors and examine the stability 
of this helical phase. We show that the helical phase typically occupies most of the magnetic field- 
temperature phase diagram. However, we also find that this phase is sometimes unstable to a 
| multiple- q phase, in which both the magnitude and the phase of the order parameter are spatially 

Q_) . varying. We find the position of this helical to multiple- q phase transition. We further examine the 

^"""^ ■ density of states and identify features unique to the helical phase. 

OO ! PACS numbers: 74.20-z, 71.18. 
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V ' T. INTRODUCTION 

Recently, a new class of heavy fermion superconductors have been discovered that break parity symmetryi 2 ^. The 
^ broken parity symmetry implies that the three materials CePtaSi, CeRhSi3, and CeIrSi3 all allow a Rashba spin-orbit 
coupling. The energy scale of this coupling is much larger than the superconducting energy scaled. This has many 
non-trivial implications on the resulting superconducting state^^ ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 . Of particular relevance 
to the work presented here, are the unusually high upper critical fields in these three Rashba superconductors. 
| These fields substantially exceed the usual Pauli limiting field. Consequently, the Zeeman interaction must play an 
important role in the the physics of the superconducting state. Some progress has been made in identifying the ground 
states of Rashba superconductors in Zeeman fields. In particular, microscopic arguments have been presented that 
indicate that a superconducting 'stripe' phase appears in two-dimensional (2D) s-wave Rashba superconductors when 
a Zeeman field is applied in the pian o 9 ' 11 . This phase resembles the Zeeman field induced phase for conventional 
( \ superconductors found theoretically by Fulde, Ferrell, Larkin, and Ovchinnikov (FFLO) 19 i 20 which has recently been 
observed in CeCoIn j 21 i 22 . In the FFLO phase, the superconducting order vanishes periodically along a line, taking 
^^O ' the simple form A(J2) = Ao cos(q • R) near the upper critical field. As the field is reduced, more Fourier components 
appear in the order parameter of this phase and, at low enough fields, a true stripe like order develops 2 ^. We 
will call this phase (and its appropriate generalization to Rashba superconductors) the multiple- q phase. However, 
phenomenological arguments for Rashba superconductors indicate that a helical phase is the stable ground state near 
\Q ' T c , once a magnetic field is applied 8 - ' 16 ' 17 . In this phase, the order parameter takes the form A (J?) = Aoe^-^. The 
gap magnitude is spatially homogeneous and therefore the helical phase exhibits quite different physical properties 
than that of the multiple- q phase. To address the physical properties of Rashba superconductors in magnetic fields, 
it is important to understand which of these phases are stabilized. 

Vortices will also play an important role in understanding the physics of Rashba superconductors in magnetic fields. 
This work focuses on the role of a Zeeman field on Rashba superconductors. Previous experience has shown that the 
physics associated with the Zeeman field persists when vortices are present. In particular, the Zeeman field induced 
FFLO phase coexists with vortices in the heavy fermion superconductor CeCoIn j 21 i 22 . Furthermore, the helical phase 
discussed above has been shown to coexist with vortices in Rashba superconductors near the upper critical fiel d 16 ' 17 . 

To understand the microscopic origin of the helical and multiple- q phases, it is useful to consider the quasiparticle 
states when inversion symmetry is broken and a Zeeman field is applied. Broken inversion symmetry allows an 
antisymmetric spin-orbit coupling: ag(k) Sf t for a quasiparticle with spin and momentum k. Note that g^ — — g ^ 
due to time reversal invariance. For a Rashba interaction, g^ = (k y , —k x , Q)/kp- The addition of the Zeeman field leads 
to the additional coupling hbH ■ Sj, and the quasiparticle energy eigenvalues become Ey, ± = £^ ± \ag^ +(J,bH\ where 
£k is th e quasiparticle energy when a = fisH = 0. We are interested in Zeeman fields that are comparable to the gap 
energy scale and therefore consider the limit a >> [IbH- For a cylindrical Fermi surface, with a Rashba interaction, 
applying a field along the x direction, we find , = kl/2m + (k y ± g/2) 2 /2m ± \ak\/kF where q = 2mfiBH/\k\. The 
key point is that the Fermi surfaces remain circular and the centers of the two Fermi surfaces are shifted along y in 
opposite directions. In this situation, one of the bands can gain condensation energy by pairing fermions through the 
new center of the appropriate Fermi surface. This leads to the helical phase in which the condensate wavefunction 
becomes A(U) = A exp(zq • R). This situation is depicted in Fig. 1. However, the gain in condensation energy 
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FIG. 1: A magnetic field directed as shown shifts the center of the two Fermi surfaces by ±q/2. The smaller dot represent the 
point (0,0) (center of Fermi surfaces without field) and the two larger dots represent the points (0, — q/2) and (0,q/2) (these 
are the centers of the new Fermi surfaces). To gain condensation energy, pairing occurs between states of k + q/2 and —k + q/2, 
leading to a gap function that has a spatial variation A(i?) = Ao exp(iq ■ R). 



of one Fermi surface is accompanied by a corresponding loss in condensation energy on the other Fermi surface, 
since the centers are shifted in opposite directions. This leads to a competition between the stability of the helical 
phase and that of the muliple-g phase (in which both ±q gap function modes appear). In this paper, we address the 
resulting phase diagram by developing the quasi-classical Eilenberger equations for a Rashba superconductor in the 
limit a >> /is-ff. We further examine the density of states of the helical phase. 



II. MICROSCOPIC FORMULATION 



We consider the following Hamiltonian: 

H ° = J2 C k,s + [a9k + ^ bH] ' ° k ^' + F(fc ' k 'K + q/2,/-k+q/2,i C -k'+q/2,l c k'+q/2A- ( l ) 

k,s,s' k.k.q 

We also set (g?) = 1 where () represents an average over the Fermi surface. We have restricted ourselves to spin-singlet 
pairing interactions since this is sufficient to capture the new physics associated with the helical phase. In the large 
a limit, a » T c , the pairing problem becomes a real two-band problem in the diagonal spinor (±) basis. In this 
basis, the pairing interaction becomes (this is after redefining the gap functions and the anomalous propagators by a 
k dependent phase factor—), 

V=\v{k,k)^\ - 1 ) . (2) 

We work within the quasi-classical approximation and define the usual Green's functions in Nambu space for each 
band & ± {R) = [ip±(R),ip±(R)] and define the imaginary time Green's function 

G ± (x 1 ,x 2 ;t 1 -t 2 ) = -(T t * ± (si,t 1 )*J : (x2,72j)), (3) 

here the operator T T arranges the field operators in ascending order of the imaginary time < r < 1/T and 
*S?(x,t) — e Tni $!(x)e~ Tn . We introduce the center-of-mass coordinate, R = (x\ +x 2 )/2 and the relative coordinate, 
r = X\ — x 2l and perform the Fourier transformation in the latter according to 

G ± (k,R;uj n ) = fdr f ' drG ± (x u x 2 ; t) e -<(*-*-« ^\ (4) 



where uj n = (2n + l)7rT is the fermionic Matsubara frequency. We define 

9± f± \ = £ 
fl -9± J ~ it 



g±(k,R,u; n ) = ( 9 i f± )=- I d£T 3 G ± (k,R,uj n ), (5) 
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where g?£ integrates out the variable perpendicular to the Fermi surface, fe is vector on the Fermi surface, and T3 is 
the z-component of the Pauli matrices acting on the particle-hole space. The standard quasi-classical approac h 23 ' 24 i 25 
results in the following Eilenberger equations for this system : 

[uj n ±ifi B gu B + v; -(V+— A)]f ± = A ± (k,R)g± (6) 



2ie 



K ± WBQ U B-v k -(V A)]f ± = A* ± (fe, R)g± 



(7) 



where f±f± +g± = l,u) n = 7rT(2n+l), fe denotes fc restricted to the Fermi surface, A is the vector potential (included 
for completeness), and B — V x A. We will neglect A and set B = H in the following. With the two band pairing 
interaction of Eq. [2J the gap equation is 



k 



(8) 



where Nj is the density of states on band j. We also express V(fe, fc ) = — V(pr(k)(f^ (fe ) where F labels the ID gap 
representation we are interested in. The gap function can be written in the product form A a (k, R) — A a (R)(pr(k) 
where Eq. [2] implies that A + (R) + A_(i?) = 0, indicating that the gap magnitudes on the two bands are equal at 
every position in space. Henceforth we set A (J?) = — A+(i2) = A_(J?). 

The Eilenberger equations can be derived from a Gibbs free-energy functional 2 ^. We will require this free energy 
to compare the energies of the different phases. Once Eqs. [5] and [7] are solved for a given functional form of A(R), 
this free energy functional becomes 



^sn = / dR 



V- 1 \A(R)\ 2 - nTj^Nj (lj(k,R, w«) 



with 



Ij(k,R,u) n ) = 



^(k)A*(R)f j (k,R,oj n )+f](k,R,u; n )A(R)>pr{k) 



gj(k,R, uj n ) 



(9) 



(10) 



III. ANALYSIS OF THE HELICAL PHASE 
A. Helical phase solution 

Eqs. [6]and[7]can be solved for A(fe, R) — A( y 5r(fc)e i ^'^, the solutions are: 



9± = , (11) 



y/<J n 2 ± + \Aip r (k)\ 2 



/± = %ffiff (12) 



and 



J n 2 ± + lA^rWI 2 

where uJ n ± = u n ± i^sg ■ H + iv ■ q. The free energy Eq. [9] is then minimized with respect to A and q to find the H-T 
phase diagram. This phase diagram depends strongly upon field orientation, dimensionality, and the relative values 
of iV_ and N + . It depends weakly upon the pairing symmetry For field parallel to z, the stable state is given by 
q = and the physics is independent of the Zeeman field 1 ' . In the following, we consider only magnetic fields applied 
perpendicular to z. When N + ^ N-, then q ^ whenever H 7^ 0. Some analytical results can be found for the upper 
critical field. For a 2D cylindrical Fermi surface, independent of pairing symmetry, the upper critical field diverges as 
T — > and Vpq = /igi x H at H C 2 for N + > iV_ (a related result has been found previously—^). We have further 
found that if both attractive spin-triplet and attractive isotropic spin-singlet pairing interactions are included, then 
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the divergence of H C 2 occurs occurs at T > for a 2D cylindrical Fermi surface. In particular, if N + = iV_, then this 
divergence occurs when T = ^/T s T p where T s (T p )is the usual T c for isotropic s-wave (p-wave) pairing when a = 0. 
For an isotropic s-wave superconductor, with a 3D spherical Fermi surface at T = 0, H C 2 = 2Hpe 1+7r \ SN \/ 2 where 
8N = (N+ - N-)/(N+ + N-), and H P = -j=^ is the usual Pauli paramagnetic field , and A = 1.76T C . Generically, 
the appearance of the helical phase enhances the upper critical field well beyond the Pauli limiting field as is observed 
in CePt 3 SiiL2£, CelrSi^, and CeRhSi^. 



B. Stability of the helical phase 

As explained in the Introduction, there are physical reasons to suggest that the helical phase is not always stable. 
To examine this possibility, we set A (A) = A q e l 1^ + A p e l P'^ + A2 q - p e l C 2 Q~P s> '^ where q is the Fourier component 
that optimizes the helical phase free energy. The other two Fourier modes represent the instability of the helical phase 
to a phase we name the multiple- g phase. With this solution for A(R), we solve the Eqs. [5] and [7] perturbatively, and 
expand the free energy in Eq. [5] to second order in these two modes. 

To carry out this procedure, we re-write the Eilenberger equations in terms of Fourier components. In particular, 
we set A(it) = Yq^^Aq, ff(^) = Ylq^ ^9g> an d f(R) = Z) 9 ' JQ ( we nave suppressed the a, k and u> n 
labels for notational simplicity). The Eilenberger equations become 

u qfq = A p9i-p ( 13 ) 
p 

U- q fl=J2 A * P 9 q + P (14) 
p 

^2i9q+p9-q + fq+pf-q] = <W (15) 

q 

where we have defined io q = u> n ± ifJ-BdL ' B + v°L ' Q (the ± refer to the two different bands). We wish to solve these 
equations as a perturbation about the helical phase Aq ^= (we use gauge invariance to choose this real), keeping 
terms in the free energy up to second order in A q and A2Q- q . To carry this out, we require f q , f^_ q ,g q ^Q, and <?q- 9 
to first order in the perturbation and go and Jq + /Iq to second order in the perturbation. We label the second 

order corrections to go, /q + f_Q as go, }q + f_g and keep the labels go, Jq, f_g for the zeroth order perturbation. 
A lengthy calculation gives 



go [(2u Q u 2Q _ q + Ag)A g - A 2 Q A* Q _ q ] 

20J Q UJ q UJ2Q-q + Ag(Wg + L02Q-q) 



fq = - n — 7%7. . w T gJ . (i8) 



goA Q [u 2Q - q A q + uj q A* Q _ q ] 
9q-Q — ~o 1 A 2 7 i \i 

2ui q LU 2 Q-qUJQ + A Z Q {Ul q + UJ2Q-q) 

7 , 7i 9og q -Q{A q + A* Q _ q ) + g g q ~ Q {A* + A 2Q - q ) - 2A Q g q -Qg Q - q - fl q f q A Q - fi 2Q+q f 2 Q- q A Q 

v go^Q + A Q f Q 

and 

^9q-Q9Q-q+ /q(/q+ f- Q ) + f-qfq + f-iQ+nfiQ-q 

90 = 7T — ■ (19) 

Using these expressions in Eq. [5] leads to a free energy of the form 

a p |A p | 2 + a 2 q- P \A2q- P \ 2 + a m A p A 2q - P + a* m A* p A* 2q _ p (20) 
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FIG. 2: Phase diagram for an isotropic s-wave superconductor with a 3D spherical Fermi surface when SN = and 5N = 0.05. 
When N+ and N_ differ slightly, the helical phase dominates the phase diagram. 



The helical phase is unstable when the above free energy becomes negative for any choice of A p or A2 P _ g . We find 
numerically that the instability occurs for p = —q. Fig. 2 shows the phase diagram found for N + = iV_ (SN = 0) 
for isotropic pairing with a 3D spherical Fermi surface. In this case, the helical phase occupies only a small region 
of the phase diagram. We find that the transition from the uniform phase into the helical phase is first order for 
T/T c < 0.39. However, for T < 0.36, we find that the helical phase is itself unstable to the multiple-g phase. Fig. 2 
also shows the phase diagram found when SN = 0.05. Contrary to SN = 0, the helical phase occupies almost the 
entire phase diagram, with the multiple-g phase existing in a region at low temperature and moderate fields. The 
region occupied by the stripe phase decreases as SN increases. Once SN > 0.25, we find that the multiple-g phase 
ceases to exist. The suppression of the multiple-g phase due to an increasing SN occurs for all pairing symmetries. 
It would nevertheless be of interest to look for the helical to multiple-g phase transition in Rashba superconductors. 
Fujimoto has pointed out that Fermi liquid corrections for heavy fermion materials leads to a large enhancement of a 
magnetoelectric effect that has the same origin as the helical phasei^ii^. In our theory, this enhancement is captured 
by an increase of SN. 



C. Density of states in the helical phase 

Having established the stability of the helical phase, it is of interest to determine some of its physical properties. 
The density of states N(u>) is an important quantity for many properties and it is given by 

N(u) = Y / N J 9t( j= M == ) (21) 

k 



3 =± Ww? - |Apr(fc)|= 



where uj± = lo ± I^BQj^ 1 H + ■ q. Here we consider in detail d x 2_ y 2 pairing symmetry on a 2D cylindrical Fermi 

surface at T/T c = 0.15 (and give analogous results for isotropic s-wave pairing). Note that d x i_ y i pairing is a possible 
pairing symmetry in CePtaSi, in which line nodes have been observe d 27 i 28 . We set SN = 0.25 and have checked that 
the helical phase is the stable ground state. For T/T c = 0.15, H C 2 = 21Hp with the field applied along the gap 
maxima. The evolution of N(lo) reveals two properties of interest that are revealed in Fig. [3J The most significant 
property is that for fields H/H C 2 >~ 0.25, N(ui) is field independent and has two contributions. The first is a normal 
component that exists on one of the bands and the second contribution is the usual d-wave density of states from the 
second band. Similar behavior exists for isotropic s-wave pairing as is shown in Fig. 3J The observation that N(u>) 
for the one of the bands corresponds to that of a normal metal is intriguing because it exists even though both bands 
have the same pairing amplitudes. Nevertheless, this result is intuitive because one band prefers VFq — (J-bz x H and 
is not frustrated while the other band prefers vpq = —\ibZ x H and therefore cannot pair the fermions that are on the 
Fermi surface. This manifests itself as an increase of N(u>). It is natural to expect that the frustrated band will have 
its gap amplitude become zero at sufficiently large fields. This possibility is not permitted by the theory considered 
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FIG. 3: Density of states in the helical phases for increasing magnetic fields for a d-wave superconductor. The high field density 
of states is field independent and corresponds to that of an ungapped band and a gapped band. 
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FIG. 4: Density of states in the helical phases for increasing magnetic fields for a s-wave superconductor. 



here because only spin-singlet pairing interactions have been included. Introducing an attractive spin-triplet pairing 
interaction will lead to a gap that vanishes on the frustrated band as the field is increased. However, even if the 
frustrated band has zero gap amplitude, the single particle density of states will remain qualitatively the same. 

The second property of interest occurs at low fields. In particular, the single peak at u> = A of the H — density of 
states splits into five peaks: one at lu = A and the others at approximately lu = |A± (IbH ±vp q\. This occurs for the 
field applied along the directions in which the gap is maximal. This property is anisotropic, for the field applied along 
the gap nodes the oj = A peak splits into four features that occur at approximately co = |A ± ± VFq/V%\- 

This property may provide a means to probe the gap symmetry in CePt 3 Si. Note that if q = 0, then an anisotropy 
still exists in N(lu). At low fields, the s-wave density of states has four features for all in-plane field orientations as is 
shown if Fig. [U 



IV. CONCLUSIONS 



We have derived the quasi-classical Eilenberger theory describing Rashba superconductors. Using this theory, we 
have examined the stability of helical and multiple- q phases for Rashba superconductors in magnetic fields. We 
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have found that the helical phase is stable over a wide range of the phase diagram. Finally, we have examined the 
single-particle density of states in the helical phase and have revealed that at large fields the it behaves qualitatively 
differently for the two spin-split bands. On one band the density of states is completely normal, while on the other 
band it remains gapped and resembles the zero-field superconducting density of states. 

We are grateful to P.A. Frigcri, N. Hayashi, and M.Sigrist for many helpful discussions. This work was supported 
by the National Science Foundation grant No. DMR-0381665. 
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